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Abstract 

We consider multiply periodic functions, sometimes called Abelian functions, defined with respect to 
the period matrices associated with classes of algebraic curves. We realise them as generalisations of 
the Weierstrafi p-function using two different approaches. These functions arise naturally as solutions 
to some of the important equations of mathematical physics and their differential equations, addition 
formulae, and applications have all been recent topics of study. 

The first approach discussed sees the functions defined as logarithmic derivatives of the cr-function, 
a modified Riemann 6'-function. We can make use of known properties of the sigma function to derive 
power series expansions and in turn the properties mentioned above. This approach has been extended 
to a wide range of non hyperelliptic and higher genus curves and an overview of recent results is given. 

The second approach defines the functions algebraically, after first modifying the curve into its equiv- 
ariant form. This approach allows the use of representation theory to derive a range of results at lower 
computational cost. We discuss the development of this theory for hyperelliptic curves and how it may 
be extended in the future. 

1 Introduction 

In this paper we work with generalisations of the Weierstrafi p-function. This is an elliptic function, and so 
periodic with respect to two independent complex periods uii,ui2'. 

p{u + uji) = p{u + UJ2) = p{u), for all u E C. (1) 

The p-function has the simplest possible pole structure for an elliptic function and satisfies many interesting 
properties. For example, the p-function can be used to parametrise an elliptic curve, 

y"^ = - g2X - gs, (2) 

where 52 and 53 are constants. It also satisfies the following well-known differential equations, 

{p'{u)Y ^ 4p(w)3 -£,2p(w) -53, (3) 
p"(u) = 6p(u)2-ig2. (4) 

Elliptic functions have been the subject of much study since their discovery and have been extensively used 
to enumerate solutions of non-linear wave equations. They occur in many physical applications; traditionally 
the arc-length of the Icmniscatc and the dynamics of spherical pendulums, [28], but also in cryptography, 
[30] . and soliton solutions to the KdV equation, [25]. 

Recent times have seen a revival of interest in the theory of their generalisations, the Abelian functions, 
which have multiple independent periods, or more accurately, are periodic with respect to a multi-dimensional 
period lattice. The lattice is usually defined in association with an underlying algebraic curve. These functions 
are also beginning to find a wide range of applications. For example, they give further solutions to the KdV 
equation along with solutions to other integrable equations from the KP-hicrarchy, (see for example [12] , [8] 
[22]). They have also been used to describe geodesic motions in certain space-time metrics, [24] . 

This introduction will continue with a discussion of the motivation to study such functions followed by a 
reminder of the Riemann-Roch theorem and its implications in this context. Then in Section [2] we discuss 
how the theory of higher genus functions can be derived by defining them via an auxiliary function, the 
(T-function. We describe how a series expansion of this function can be derived and the results which follow. 
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In Section |3] we proceed to give an alternative definition of the functions, developed in recent years. This 
algebraic definition and the following approach is termed equivariant and brings out an s^2(C) structure in 
the equations which can be exploited to reduce the computation involved in working with the functions. 
Finally in Section |4] we discuss the possible future directions of this research. 



1.1 Motivation 

As discussed above, the Weierstrafi p-function and its generalisations have strong connections with integrable 
PDEs. The KdV equation has a translational symmetry reduction described by the Weierstrafi p-function 
and hence is associated with a genus one curve. More generally the KdV with other specific, differential 
constraints is associated with generalised p-functions and genus g hyperelliptic curves. For details see the 
example below. Of course, such solutions also arise via the Integrable Scattering Transform and are expressed 
in terms of ^-functions. So the interest for Integrable Systems is, in part, to obtain equations satisfied by 
p-functions for general plane algebraic curves. On a wider plain there is interest in application to problems 
in GR, Billiards, Statistical mechanics and elsewhere. Most fundamentally we are concerned here with the 
interplay between function theory and the geometry of curves. 

Example: A genus two curve 

Consider for x,y G C, the genus two curve 

y2 = Agx*' + 6A5X^ + IbXiX^ + 20X3X^ + ISAaa;^ + 6Xix + Aq. 

In this case there are generalised p-functions, denoted pij{ui,U2) for i,j = 1,2. They are functions of g = 2 
variables and are multiply periodic. We can denote the derivatives of the functions using additional indices, 
so for example 

d 

For simplicity we shorten the differentiation notation so that 

dui 

for example. It is the case that pij — pji and so we have three generalised p-functions. We also have 

dipjk = djpik- 

and so generally the order of the indices is not relevant. We adopt the notation of writing the indices in 
ascending numerical order. 

These functions can be shown to satisfy the following set of partial differential equations. 

-3P2222 + 2p22 = A4P22 - 2A5P12 + Aepii + A2A6 - 4A3A5 + 3A4 
-5P1222 + 2pi2p22 = A3P22 - 2A4P12 + A5P11 -I- ^(AiAe - 3A2A5 + 2A3A4) 

-^Pll22 + f P11P22 + |p?2 = •^2p22 - 2A3P12 + A4P11 + i(AoA6 - 9A2A4 + 8A^) 

-3P1112 + 2piipi2 = A1P22 - 2A2P12 + A3P11 + i(AoA5 - 3A1A4 + 2A2A3) 



-|piiii + 2pfi = A0P22 - 2A1P12 -I- A2P11 -t- A0A4 - 4A1A3 -t- 3A 



Consider the first equation. By a rational transformation we may, without loss of generality, set Ag = 0, 
A4 = 0, A5 = i to give 



1 o 2 1 2^ 

-^P2222 + 2p22 = -gPl2 - ^A3. 



We differentiate the identity with respect to U2 to find 

1 , 1 

— 2P22222 +4p22p222 — — 2P122, 

and making the change of variables di = dt, 82 = dx P22 = U{x, t) this is just 

Ut - U^.^.^, + \2UUx = 0. 
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Thus wc have the KdV-equation, along with four additional differential identities. 

We can write these five partial differential equations in a parametrised form. Thus the relations on 4-index 
objects are 

- Id^p + 2{d^p)^ = HoB^p - 2HiBdp + H^ff^p + G (5) 

where d = di + vd2 , B ~ 82, the Hi arc functions of Ai and G is a function of A^ and i'. Another set of identities 
in the function are those quadratic in the 3-index objects. These may be written in the parametrised form 

d^pd^p - 6id^p)iddp){d^p) = H'iB^p, Bdp, d^p) 

where H' is a polynomial of degree two in the second derivatives of p with coefficients in A^ . 

Equation ([S]) is a reduction of the Boussinesq equation and so we can write down a Lax pair for the genus 
two p-function equations: 

M = d^ + iC'a^ + ^,{eyd - 3{d^p)d - lid^p) - Cid^p) + ICidBp) 
where = .g(i.) and C = f , SI- 
1.2 Riemann-Roch Theorem 

Let X be a compact, non-singular, complex algebraic curve (Riemann Surface) and D a divisor on X: 

D = niPi + n2P2 + . . . + UrPr, Pi eX,nie Z. 

The Riemann-Roch theorem relates dimensions of spaces of functions and of holomorphic differentials asso- 
ciated with a divisor to the degree of the divisor and the genus of the curve: 

dim H*^ - dim H^ = 1 - g + dcgD. 

Here i/^ is the C- vector space of meromorphic functions having poles of order at most rii at each Pi, Hj^ 
the C- vector space of holomorphic differentials with zeros of order at least rii at each Pi. 
More useful for our immediate purpose is the formulation, 

(dim iJ^+P - dimi?^) + (dimiJ^ - dimi?|j+p) = 1. 

Thus it is that in adding a point P to the divisor D, either H'^ augments by +1 and H^ is unaltered or 
H^ augments by —1 and H'^ remains unaltered. 

At a generic point, P, on the curve, consider the family of divisors nP for n e N U {0}. The following 
pattern of dimensions obtains. 

D P ... gP {g + l)P {g + 2)P ... 
dimi/^ 1 1 ... 1 2 3 

dimHl g g-1 ... 

Note that there are g 'gaps' where the addition of an extra P results in no new functions on X. Each gap 
is associated with the loss of a holomorphic differential. The pattern is, however, broken at special choices of 
divisor, more specifically special choices of P in the sequence above. At such Weierstrafi points the interval 
of g gaps is broken as illustrated in the examples below. 

Example: Genus two hyperelliptic case 

There are Weierstrafi gaps at P and 3P: 



D 





P 


2P 


3P 


4P 


5P 


6P 


IP 


8P 


9P 


lOP 


dimi/^ 


1 


1 


2 


2 


3 


4 


5 


6 


7 


8 


9 
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The curve is a relation on i?iop' 



The notation [x"^] denotes a polynomial of degree n in x. A relation arises at lOP because we can create 
elements in H'^Qp\Hgp in two ways when only one extra dimension is available. The diagram also indicates 
the form of holomorphic differential (depending on f{x, y)) lost at each gap. 

Example: Genus three non-hyperelliptic case 

There are three Weierstrafi gaps at P, 2P and 5P. 
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The curve is a relation on H^2P'- 



f{x,y)^y^ + [x]y^ + [x^]y+[x^]^0. 



fix,y) 



Transformations 

In each of the examples above x and y are coordinates on X, i.e. maps X —i' F^. We shall be interested in 
coordinate transformations 

x ^ X = 'ip o X for : — !• . 

Note that under such transformations the coordinates x and x have poles at distinct divisors on X: Weierstrafi 
points for x map to more general special divisors for x and give different models of X: rationally equivalent, 
and usually singular, curves in P^. For example, the following three curves are all rationally equivalent. 

y^ + [x^][xf^O ^ f + [x^][S:]^ = ^ [x]f + [x^]=0 

It so happens that this example is of genus three. Every genus three X has a nonsingular embedding in P^ 
and the third equation above is it, j26j . 



2 Working with Kleinian functions 

Weierstrafi introduced an auxiliary function, cr(u), in his theory which satisfies 

cP 

The cr-function plays a crucial role in both the generalisation and in applications of the theory. It satisfies 
the following two term addition formula, 

a{u + v)a{u - v) , ^ , ^ 

The Weierstrafi a-function is usually defined as an infinite product over the periods, but it can also be 
expressed using the first Jacobi ^-function, multiplied by a constant and exponential factor. It is this 
definition that generalises naturally to give higher genus cr-functions as multivariate functions defined using 
the Riemann 0-function, as presented in Definition 12.31 Once the generalised cr-function is defined we may 
then consider generalised p- functions in analogy to equation ([6]), as presented in Definition 12.51 

This approach was pioneered by Klein and Baker as described in Baker's classic texts [5] and [7]. Hence 
these generalised functions are sometimes called Kleinian. The periodicity conditions of all the higher genus 
functions are defined using the period lattice of an underlying algebraic curve. In this section we work with 
the following class of curves, while in Section |3] we consider an equivariant class. 
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Definition 2.1. For two coprime integers {n,s) with s > n we define an {n, s)-curve, denoted C, as an 
algebraic curve defined by f{x, y) ~ 0, where 

f{x,y) = 2/" + p2ix)y''-^ + ■ ■ ■ + Pn-iix)y ~ p„{x). (8) 

Here x, y are complex variables and Pj{x) are polynomials in x of degree (at most) ljs/n\. We define a 
simple subclass of the curves by setting Pj{x) — for < j < n — 1. These curves are defined by 

fix, y) = - (x'' + Xs-ix'-^ + • ■ • + Alt + Aq) (9) 

and are called cyclic {n, s) -curves. We denote the curve constants by Xj for the cyclic curves and by fij 
for the general (n, s)-curves. Note that in the literature the word 'cyclic' is sometimes replaced by 'strictly' 
or 'purely n-gonal'. 

Restricting to the cyclic classes results in much easier computation, but it not usually necessary for 
theoretical reasons. However, the cyclic curves do possess extra symmetry which manifests itself in the 
associated functions satisfying a wider set of addition formulae. 

If the curve C is a non-singular model of X then the genus is given by g — ^{n — l)(s — 1) and the 
associated functions, to be defined shortly, will be multivariate with g variables, u = (ui,...,Ug). As an 
example, the elliptic curve in equation ([2]) is a (2,3)-curve and the associated Weierstrafi a and p-functions 
depend upon a single complex variable u. 



2.1 Defining the functions 

We first describe how the period lattice associated to the curve may be constructed. We start by choosing 
a basis for the space of differential forms of the first kind; the differential 1-forms which are holomorphic on 
the curve, C. There is a standard procedure to construct this basis for an (n, s)-curvc, (see for example [29]). 
In general the basis is given by hidx/ fy, i = 1 ... 5 where the hi are monomials in (x, y) whose structure may 
be predicted by the Weierstrafi gap sequence, although other normalisations are sometimes used. 

We next choose a symplectic basis in Hi (C, Z) of cycles (closed paths) upon the compact Riemann surface 
defined by C. We denote these by {ai, . . . a^, . . . We ensure the cycles have intersection numbers 

Ui ■ aj = 0, f3i ■ =0, tti • ft = 5ij = 

The choice of these cycles is not unique, but the functions will be independent of the choice. 

We introduce dr as a basis of differentials of the second kind. These are meromorphic differentials on C 
which have their only pole at 00. This basis is usually derived alongside the fundamental differential of the 
second kind. Rather than repeat the full details here we refer the reader to [5] for the general theory and 
which gives a detailed example construction for the (3,5)-curve. 

We can now define the standard period matrices associated to the curve as 




We define the period lattice A formed from uj' ,uj" by 

Note the comparison with equation ([T]) and that the period matrices play the role of the scalar periods in 
the elliptic case. The functions we treat are defined upon with coordinates usually expressed as 

U = (Ul, . . . ,Ug). 

The period lattice A is a lattice in the space C^. The Jacobian variety of C is presented by C^/A, and is 
denoted by J. We define k as the modulo A map. 



/ 1 iH = j, 
I iii^j. 
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For k = 1, 2, . . . define 2l{, the Abel map from the fc-th symmetric product Sym'^(C) of C to J by 
2(( : Sym''^(C) ^ J 

(Pi, . . . , Pfc) ( f du+---+ f du\ (mod A), 

y OO OO / 

where the Pi are again points upon C. Denote the image of the fc-th Abel map by W^''^ and define the k-th 
standard theta subset (often referred to as the fc-th strata) by 

where [—1] means that 

[-l]{Ul,...,Ug) = {-Ul,...,-Ug). 

We are considering functions that are periodic with respect to the lattice A. 

Definition 2.2. Let 9K(m) be a meromorphic function of u S C^. Then 9Jt is a standard Abelian function 
associated with C if it has poles only along k^^[Q^^~^^) and satisfies, for all £ G A, 

mi{u + i) = m{u). (10) 

The generalisations of the Weierstrafi p-functions with which we deal will satisfy equation ([TU| . The 
generalisation of the ct- function we define below will be quasi-periodic. Let 8 ~ u' 5' 4- oj" 8" be the Riemann 
constant with base point oo. Then [8] is the theta characteristic representing the Riemann constant for 
the curve C with respect to the base point cxd and generators {a^, jSj} of i?i(C, Z). (See for example [T2] 
pp23-24.) 

Definition 2.3. The Kleinian cr-function associated to a general (n, s)-curve is defined using a multi- 
variate 9 -function with characteristic 8 as 

a{u) = ceiip{^wf{Lo')-^u'^)-9[S]{{Lu')-'^u'^\{Lu')-^Lu"). 



cexp (^^ur]'{uj') ^m^) ^ exp 



27ri{i(m + 8'f{uj')-^uj"{m + 8') + (m + S')^ {{uj'y^u'^ + 6")} 



The constant c is dependent upon the curve parameters and the basis of cycles and is fixed later, following 
Lemma[ 



We now summarise the key properties of the cr-function. See [12] or [29] for the construction of the 
cr-function to satisfy these properties. For any point it 6 we denote by u' and u" the vectors in such 
that 

u = uj'u' +uj"u". 

Therefore a point £ G A is written as 

£ = Lu'£' + Lu"£" G A, G Z". 

For u,v G and £ <= A, define L{u, v) and x(-^) as follows: 

L{u,v) = u^{ri'v' + ,fv"), 

Xi£) = exp[2ni{i£r6"^i£"fS' + m^n]- 
Lemma 2.4. Consider the a-function associated to an {n, s)-curve. 

• It is an entire function on . 

• It has zeros of order one along the set k~^(0[s~^1). Further, we have cr{u) ^ outside the set. 

• For all u G <C^ ,£ G A the function has the quasi-periodicity property: 



a{u + £)= x{i) exp 



<7{U). 
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• It has definite parity given by 

(J{-U) = 

Proof. The function is clearly entire from the definition, while the zeros and the quasi-periodicity are classical 
results, (see [5]), that are fundamental to the definition of the function. They both follow from the properties 
of the multivariate 0- function. The parity property is given by Proposition 4(iv) in |29| . 

□ 

We can now finally define p- functions using an analogy of equation ((5]). Since there is more than one 
variable we need to be clear which we differentiate with respect to. We define multiple p-functions and use 
the following index notation. 

Definition 2.5. Define m-index Kleinian p-functions for m > 2 by 

where ii < • • • < im G {1, • • ■ , ,9}- 

The m-index p-functions are meromorphic with poles of order m when <j{u) = 0. We can check that they 
satisfy equation (jlOp and hence they are Abelian. The TO-index p-functions have definite parity with respect 
to the change of variables u — > [— This is independent of the underlying curve, with the functions odd 
if m is odd and even if m is even. Note that the ordering of the indices is irrelevant and so for simplicity we 
always order in ascending value. 

2.2 Weights and expansions 

Klein considered generalisations of the p-functions associated with certain hyperelliptic curves. These can be 
thought of as the (n, s)-curves with n ~ 2, (when s = 3 the hyperelliptic curves reduce to elliptic curves). 
The (n, s)-curves where n = 3 are labeled trigonal curves and those with n = 4 are labeled tetragonal curves. 
The original work of Klein and Baker motivated the general definitions above, which were first developed 
in [12] and [18]. Many results for hyperelliptic curves were presented in [Ij with these methods applied to 
trigonal curves by the same authors in [T3]. The tetragonal curves were first considered in [^ . 

The methods of Klein become harder to implement when n — 3 and harder still when n = 4. The 
difficulties are not just computational but theoretical. This led to the development of a related approach 
making use of the series expansion of the cr-function. This approach was applied to the two canonical trigonal 
cases, first in [17] and [8] and then again after further development in [19] and [21]. It was essential to the 
development of the tetragonal theory in [22] and the theory of higher genus trigonal curves in [20] . In this 
section we discuss how such a series expansion may be derived and used. 

Tiie Sato weights 

For a given (n, s)-curve we can define a set of weights, denoted by wt and often referred to as the Sato weights. 
We start by setting wt(a;) = — n, wt{y) = —s and then choose the weights of the curve parameters to be such 
that the curve equation is homogeneous. We see that for cyclic curves this imposes wt(Aj) = — n(s — j) while 
for the non-cyclic curves we usually label the fXj with their weight. For example, the general (2,5)-curve is 
given by 

+ y[l^ix^ + IJ'S.x + Ms) = a;^ -I- y.2X^ + IJ-4X^ + l^&x^ + IJ^sx + Mio, 
and the cyclic restriction by 

y'^ =x^ + X4x'^ + XiX^ + \2X^ + \xx + Ao. 

In both cases the weight of x and y are respectively 2 and 5. Both equations are homogeneous of weight —10 
and so we see that the coefficients y^j are labeled with the negative of their weight while the weights of the 
Ao, Ai, A2, A3, A4 are given by 10, 8, 6, 4, 2 respectively. 

The Abel map 2li gives an embedding of the curve C upon which we can define ^ = as the local 
parameter at the origin, 2li(oo). We can then express y and the basis of differentials using ^ and integrate 
to give series expansions for u. We can check the weights of u from these expansions and see that they are 
prescribed by the Weierstrafi gap sequence. For example, in the (2,5)-casc u = (mi,W2) has weight (2, 1). 



7 



By considering Definition 12.51 we see that the weight of the p- functions is the negative of the sum of the 
weights of the variables indicated by the indices, irrespective of what the weight of cr(tt) is. We note that 
curves of the same genus wih, notationally, have the same p-functions, but may exhibit different behavior 
as indicated by the different weights of the variables and functions. We will discuss the weight of the a- 
function below. All other functions discussed are constructed from a or p-functions and their weights follow 
accordingly. We can show that all the equations in the theory are homogeneous in these weights, with a more 
detailed discussion of this available, for example, in [52]. 

The sigma-function expansion 

We can construct a series expansion of the cr-function about the origin, as described below. 
Lemma 2.6. The Taylor series expansion of a{u) about the origin may be written as 

QC 

a{u) = K ■ 5W„,,(m) + Ck{u). 

k=0 

Here K is a constant, SWn,s the Schur-Weierstrafl polynomial generated by {n,s) and each Ck a finite, poly- 
nomial composed of products of monomials in u of weight k multiplied by monomials in the curve parameters 
of weight ~{wt{a) — k). 

Proof. We refer the reader to [53] for a proof of the relationship between the a-function and the Schur- 
Weierstrafi polynomials and note that this was first discussed in [T3]. We see that the remainder of the 
expansion must depend on the curve parameters and split it up into the different Ck according to the weight 
split on terms between u and either the pLi or the Ai. We can see that each Ck is finite since the number of 
possible terms with the prescribed weight properties is finite. In fact, by considering the possible weights of 
the curve coefficients we see that the index k in the sum will actually increase in multiples of n. 

□ 

The Schur-Weierstrafi polynomials are Schur polynomials generated by a Weierstrafi partition, derived in 
turn from the Weierstrafi gap sequence for (n, s). See [29] and [13] for more details on these polynomials. 
This connection with the Schur-Weierstrafi polynomials allows us to determine the weight of the cr-function 
as 

wt((T) = (l/24)(7i2 - l)(s2 _ 1). 

In Definition 12.31 we fix c to be the value that makes the constant K = \ m the above lemma. Some other 
authors working in this area may use a different constant and in general these choices are not equivalent. 
However, the constant can be seen to cancel in the definition of all Abelian functions, leaving results between 
the functions independent of c. Note that this choice of c ensures that the Kleinian cr-function matches the 
Weierstrafi cr-fimction when the (n, s)-curve is chosen to be the classic elliptic curve. We note that there are 
other definitions of the cr-function and a definition more in line with the equivariant approach discussed in 
Section |3] is currently a topic of research. 

The expansion can be constructed by considering each Ck in turn, identifying the possible terms, forming a 
series with unidentified coefficients and then determining the coefficients by ensuring the expansions satisfied 
known properties of the cr-function. For example, ensuring it vanishes on the strata or insuring that the 
known identities between the p-functions are satisfied. 

Large expansions of this type were first introduced in jllj . in which used the generalised cr-function to 
construct explicit reductions of the Benney equations, (see also [9], [10] and [23]). Since then they have been 
an integral tool in the investigation of Abelian functions. Recently computational techniques based on the 
weight properties have been used to derive much larger expansions and we refer the reader to [22] and [20] 
for a more detailed discussion of the constructions. 

We note that such expansions are possible for the general [n, s)-curves, but that the calculations involved 
are far simpler for the cyclic cases. However, even in these cases they grow quickly in size and become even 
larger as the genus grows. For example, in |22| the construction of the cr-expansion for the cyclic (4,5)-curve 
was described in detail. Here the cr-function had weight 15 and the expansion was calculated up to C59 which 
had 81,832 non-zero terms, (from 120,964 possible terms with the correct weight structure). 

The number of Cfc required will depend on how the expansion is to be used. If used to identify relations as 
discussed in Section 12.41 then it will need to contain information on the monomials in curve coefficients that 
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Space 


Dimension 


Basis 




r(o) 


1 


{1} 




r(i) 


1 


{1} 




r(2) 


4 


il, ic3n , (£3i9, Dp?! 




r(3) 


9 


{1, pll, pl2, p22, pill, Pll2, Pl22, P222, S} 




r(4) 


16 


{1, . . . ,S, pun, . . . , p2222,3iS,52S} 




r(TO) 


2 

TO 


{basis for r(m - 1)} U {{pii...j„ }, {Si^ ...9 





Table 1: Table of bases for Abelian functions associated with a genus two curve. 



may be present in the identity in question. This can be determined by the weight and type of the desired 
identity. 

2.3 Bases of Abelian functions 

We can classify the Abelian functions according to their pole structure. We denote by r(m) the vector space 
of Abelian functions defined upon J which have poles of order at most to, occurring only on the 0-divisor; 
the strata Olf^^I where the 9 and cr- functions have their zeros and the Abelian functions their poles. 

A key problem is the generation of bases for these vector spaces. Note that the dimension of the space 
r(77i) is TO^ by the Riemann-Roch theorem for Abelian varieties, (see for example [23)- The first step in 
constructing a basis is to include the entries in the subsequent basis for r(TO — 1). Subsequently only functions 
with poles of order exactly to need to be sought. The Kleinian p- functions are natural candidates and indeed 
are sufficient to solve the problem in the elliptic case. Here g = 1 and so only one new function is required 
at each stage, which can be filled by the repeated derivatives of the WeierstraB p-functions. 

However, if 17 > 1 then new classes of functions are required to complete the bases. We discuss the genus 
two case as it illuminates the connections between the two approaches to p-functions discussed in this paper. 
Table [T] shows how bases in this case may be constructed for arbitrary m. As in the genus one case the basis 
for r(2) may be formed using a constant and the 2-index p-functions. However, when considering r(3) we 
find that even after including the 3-index p-functions we still need an additional function. The function 

5 = P11P22 - P?2 

is usually taken to fill this hole. The individual terms in S have poles of order 4, but when taken together 
they cancel to leave poles of order 3 as can easily be checked using Definition 12.51 We then proceed to 
consider r(4) and find that two new functions are needed after the inclusion of the 4-index p-functions. The 
two derivatives of S can play this role and we use the notation 




for simplicity, (and similarly for other functions). Naturally the derivatives of the S have poles of order four 
but we must also ensure that they are linearly independent of the other functions. This can be checked 
trivially by noting that they are of a different weight to both each other and all the other elements of the 
basis. 

It is simple to check that when considering higher values of to we find that the preceding basis and its 
unique derivatives always give the required number of functions for a new basis, and that they are all of 
unique weight and so linearly independent. Hence the general basis is as described in Table [l] where {•} is 
indicating all functions of this form. 

The genus 1 and 2 cases, in which we get to a stage where new bases are calculated from the old ones 
and their derivatives, are special. They fall into the class where the theta divisor is non-singular and the 
2?- module structure of such cases is discussed in |16j . However, subsequent (n, s)-curves are not covered by 
this case and so new methods must be used to derive bases here. Two approaches have been developed to find 
additional functions. (Note that in each case we still need to check the linear independence of the functions, 
which can be done using the cr-cxpansion.) 
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The first is to match poles in algebraic combinations of p-functions so that they cancel, analogous to 
5 G r(3). For example, consider the following function, constructed of products of 2 and 3-index p-functions. 

Bijklm ~ PijPklm + ^{pjkPilm + PjlPikm + PjmPikl — '^PklPijm — '^PkmPijl — "^PlmPijk) ■ 

Each term here has poles of order five, but overall both these and the poles of order four cancel (as can easily 
be checked using Definition I2.5|) and so the function belongs to r(3). This particular function can be used to 
construct classes of differential equations bilinear in the 2 and 3-index p-functions, discussed in more detail 
in |21| . This paper also constructs other classes of such combinations in a systematic way and uses them in 
an example case to find new addition formulae and differential equations. (See also |19j for more results that 
followed from this approach.) 

The second approach is to define Abelian functions by the application of differential operators on a- 
functions. For example, define Hirota's bilinear operator as Vi = d/dui — d/dvi. Then define the m-index 
Q-functions (for m even) as 

ii < ••• < in e {1, • • • ,.9}- 

These are all Abelian functions with poles of order two and so can be used to solve the basis problem in 
general for r(2). 

Both these approaches construct new functions that have a given pole structure in general, that is without 
reference to a specific underlying curve. However, there are some interesting examples of functions which 
have poles of lower order on certain curves. Consider the following A-function, originally introduced by 
Baker in [6]. 

A = pllp33 - Pl2p23 - Pl3 + Pl3p22- 

The terms in this expression have poles of order four and overall the expression will have poles of order 
three. (This may be checked using Definition 12. 5|) . However, in the case of the (2,7)-curve, these can be 
shown to cancel to leave poles of order two. The (2,7)-case has g = 3 and so 2*^ = 8 functions are required 
to form a basis for r(2). Such a basis may be formed from a constant, the six 2-index p-functions and the 
A-function. Similar functions (linear combinations of quadratics in 2-index p-functions) were constructed in 
[l9] to complete the basis of r(2) in the (2,9)-case and it appears that such functions are a feature of the 
hyperelliptic cases. Their appearance is more natural when considering the equivariant approach as discussed 
in Section [31 

2.4 Deriving differential equations and addition formulae 

The cr-expansion and bases of functions discussed above may be used to derive various identities and properties 
of the functions. For every (n, s)-curve there exist sets of differential equations that generalise the elliptic 
identities (jS]) and (j?]). However, the exact form of these generalisations is not clear. In the hyperelliptic cases 
they have an analogous form, so that the generalisation of 

p"iu) = 6piu)^ - \g2, 

is a set of equations expressing each 4-index p-function as a quadratic equation in 2-index p-functions. 
Similarly the generalisation of 

(p'(m)) =4p(u)^^-.g2p(u)-ff3, 

is a set of equations expressing products of 3-index p-function as a cubic equation in 2-index p-functions. 
However, in non-hyperelliptic cases it is not possible to find complete sets, or more accurately, we need to 
use all the functions in the basis of r(2) to construct the equations, rather than just the 2-index p-functions. 
These and other sets of differential equations between the p-functions are discussed in detail in [21]. It is 
expected that the precise structure of these equations will become apparent following development of the 
equivariant approach. 

In low genus cases such differential equations follow from an algebraic definition of the p-functions as 
described in |12j . However, for higher genus cases the cr-expansion is necessary to find an exhaustive set. 
Individual relations can be calculated by forming the relations with unidentified coefficients, based on the 
weight and structure of the relations. The coefficients may then be identified by substituting in the a- 
expansion, setting the coefficients to zero and solving the overdetermined resulting linear system. The 
procedure to do this is discussed in more detail in [52] and [5D]. 
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A similar use of the cr-cxpansion is in the construction of addition formulae. The Weierstrafi formula in 
equation ([7]) has been generalised in many cases and papers. The simplest generalisation is to a formula of 
the form 

(t(u + v)a(u — v) v-^ . , , 

where A^, Bi g r(2) and the c,; are constants. 

This structure follows from linear algebra after checking that the left hand side is both Abelian in u and 
V. The coefficients q can be identified by forming the right hand side from elements of r(2) and arbitrary 
coefficients and then substituting in the cr-expansion. There are many simplifications to be made to this cal- 
culation by considering the weight, parity and symmetry properties of the left hand side and constructing the 
right hand side accordingly. Further computational simplifications can be made by expanding the resulting 
products of series in a way that takes into account the homogeneity of the weight structure. 

There are a second class of addition formulae present only when the underlying (n, s)-curvc is cyclic, 

y" = x'* + A^-i-T^^i + ... + Aia; + Ao. 
It is associated with their extra family of automorphisms: 

[C-']:(x,2/)^(x,C^y), where C" = 1, 
with the following function becoming Abelian in the n variables m^, i = 1 . . . n. 

A A (t((mW)") 

It can hence be written as a sum of products of n functions from r(n) in the variables, mW,^ = 
respectively. 

Such formulae can be computationally difficult to compute. Simplified versions may be found, setting one 
or more of the variables to zero. We can also consider reduced curves which have further automorphisms 
and hence extra addition formulae. For example if we set all the curve coefficients except Aq to zero then 
there will be automorphisms 

\(f \ : {x,y) I— (r]-'x,y), where = 1. 

The first example of these automorphism addition formulae was given in [17j with further examples 
recently published in [19] and [21]. 



3 The equivariant approach 
3.1 Equivariant curves and functions 

The motivations behind the equivariant theory are to shortcut some of the calculational intricacies of the 
cr-expansion and to keep a firm grip of the underlying geometry of the curve. This approach was the topic 
of [1] where it was developed for an equivariant hyperelliptic curve of genus which is given by 




We start by choosing simple PSL{2,C) coordinate transformations on the curve coordinates 

ax + (3 y 

X X — r, y ^ y — 7 ■ 

^x -\- (7a; + OjP 

Then Hq, (as defined in Section rO)) is a g-dimcnsional S'L(2,C) module which may be reducible. It should 
be noted that this choice of transformation does not exhaust all possible and possibly useful transformations. 
They provide a convenient and manageable starting point. 



11 



We define functions pij associated with a hypcrcUiptic curve of genus g algebraically as those satisfying 

1 t-i 3-t _ VVk - F{x,Xk) 

Here the quantity / = j/j/^ — F{x, Xk) is an 'cquivariant polar form' of the curve. For example, if the genus 
is three then 

F{x, Xk) = \s{xxk)'^ + iXrixxkfix + Xk) + XeixxkfiSx'^ + 8xxk + 5x1) + i\5{xxk){x^ + Qx^xu + Qxxl + xl) 
+ 4:X4{x^ + Wx^Xk + SGx'^xl + Wxxl + x^) + 4:X3{x^ + Qx^Xk + &xxi + xl) 
+ X2i3x'^ + 8xxk + 3x1) + 4Xi{x + Xk) + Aq. 

This definition guarantees the equivariance of the algebraic definition at the expense of some small changes 
(addition of constants) to the classical definition of the pij . 

We denote the standard, finite dimensional irreducible SL{2, C)-modulc of dimension d as d, and use 
similar notation for other modules. In the hyperelliptic (2, 2g + 1) cases (of genus g), Hq is a g and the pij 
decompose as 

g g = (2g - 1) © (2g - 5) © . . . 

For the trigonal (3, 3p + 1) cases (of genus Bp), Hq is a p © 2p and the p.y decompose as 

p p = (2p - 1) © (2p - 5) © . . . 
p 2p = (3p - 1) © (3p - 3) © (3p - 5) . . . 
2p 2p = (4p - 1) © (4p - 5) © . . . 



3.2 Identities 

Products of p-functions are to be thought of as symmetric tensor products of finite dimensional modules. 
Then any identity must belong to some finite dimensional, irreducible module of identities classified by some 
highest weight element. The strategy then is to: 

• start from the Kleinian, algebraic definition of pij] 

• undertake a (classical) singularity expansion, analogous to the classical elliptic methods, to obtain some 
highest weight identity; 

• generate the module of identities associated with this highest weight identity; 

• execute a further singularity expansion modulo this module to obtain further identities. 

The tools are provided by the representation theory: the use of the Lie algebra sl{2, C) to locate highest 
weight elements and the Casimir element of the enveloping algebra to identify dimensions. 



Example: Explicit Analogue to the Weierstrafi cubic in the genus three hyperelliptic case 

Define the following matrices: 



P 






Pii 

Pl2 
Pl3 





-2pii 
-pi2 

p22 - 2pi3 
P23 



A 



Pii 

-Pl2 
2pl3 - 2p22 
-P23 
P33 



Pl2 
p22 - 2pi3 
-P23 
-2p33 





-p333 







P333 

~P233 Pl33 
P223 — Pl33 — Pl23 

~P222 + 2pi23 pl22 " Pll3 



Pl3 
P23 
P33 




P233 
-Pl33 



Pll3 
-Pll2 



H 



-p223 + pl33 
Pl23 
-Pll3 


pill 



Ao 4Ai 6A2 4A3 A4 

4Ai I6A2 24A3 I6A4 4A5 

6A2 24A3 36A4 24A5 6A6 

I6A4 24A5 16A6 4A7 

4A5 6Xe 4A7 As 



4A3 
A4 



P222 — 2pi23 
— pV22 + P1I3 
P1I2 
-pill 
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Then all identities of degree two in 3-index p-functions and degree three in 2-indcx p-functions are contained 
within, 

H-2P F 



(FAk)(l'^Ak') + -dct 



1' 
k' 








= 



(11) 



where the k, I', k' are vectors of dimension 5 with arbitrary parameters, [2]- 



For comparison, the genus one equation is written below in a corresponding fashion. (Note that here 



Pill = p')- 



- dot 
4 



A2 - 2pii 
2A3 
A4 



= 0. 



Ao 2Ai 

2Ai 4A2 + 4pii 
A2 - 2pii 2A3 

Regarding the genus 3 case we make some remarks. Firstly, the matrix A is of rank two. As a 5 x 5 
anti-symmetric matrix this implies 'Pliicker-like' relations on the pijk- This correspondence between the 
entries of A and Pliicker coordinates on the Grassmanian G(2,5) can be made precise and as an expression 
of Koszul duality which, in turn, suggests a natural conjecture for higher genus identities in the hyperelliptic 
cases. 

Secondly, equivariant non-hyperelliptic cases are still under construction. The known formulae are cer- 
tainly equivariant at leading order in the pi-^i2...i„ but not being derived from fully equivariant models of X 
the devil is in the lower order detail. 



4 Combining the approaches 

We aim to fuse the two approaches. We already know how to map from known hyperelliptic equivariant results 
to the standard results for (n, s)-curves. For example, consider the equivariant theory for the hyperelliptic 
genus 3 curve discussed in Section [X^ The set of identities encoded in the matrix equation ([TT|) is equivalent 
to the set derived and presented for the cyclic (2,7)-curve in Theorem 4.1 of [19]. To map from the former 
to the latter we need to shift the 2-index p-functions as follows: 

Pl,l Pl,l + 3A2, P1.2 P1.2 + 2A3, pl,3 Pl,3 + iA4, 
P2,2 ^ P2,2 + 9A4, p2,3 ^ P2.3 + 2A5, p3_3 p3_3 -t- SAg. 

Then we need to change the scaling of the curve coefhcients: 

Ao ^ 4Ao, Ai M> iAi, A2 t-!> iA2, A3 i-> j^Xa, A4 1-^ i^A4, 
A5 j^As, Ae ^ iAe, A7 i, ^ 0. 

Since the 2-index p-function only shift by a constant we find that the higher index p-functions are identical 
in both cases. 

The equivariant technology can help make sense of the structure of the pole bases of the functions discussed 
in Section [231 Consider again Table[T]which summarised the bases in the genus two case. The starting point 
was the three 2-index p-functions, which formed a 3-dimensional representation in the equivariant theory. 

3 = {pii, Pl2, p22}. 

If we consider tensoring this with itself then we find 

3 3 = 5 © 1, 

so we have a 5-dimensional representation and a 1-dimensional representation of functions quadratic in the 
2-index p-functions. We have 

1 = S 

= P11P22 - pL e r(3). 

The other representation is 

5 = {pIi, 4piipi2, 4pi2 + 2piip22, 4pi2p22, P22} 
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Space 


Dimension 


Basis 


r(o) 


1 


{1} 


r(i) 


1 


{1} 


r(2) 


4 


{1,3} 


r(3) 


9 


{l,3,a3,7ri[3©3] =S} 


r(4) 


16 


{l,...,E,d'3,dE} 


r(m) 


m 


{basis for r(m - 1)} U {d"'+^3,d"'-^E} 



Table 2: Table of bases for Abelian functions associated with a genus two curve. 



which are all linearly independent functions with poles of order 4 and so can replace the 4-index p-functions 
in r(4). So tensoring up the representations allows us to move through the bases and see how they break up 
into representations themselves, often in ways not apparent from the notation used. 

We can also gain insight and move between the representations by considering the representation of 
operators 

d = {di,d2}. 

Since this is two dimensional we will have 

c?(g)n = n — l©n + l 

when it is tensorcd with any representation of functions. For example, tensoring it with the 3 from above 
gives 

903 = 2®4 

where 4 is the representation of 3-indcx p-functions and 2 is empty due to the intcgrability properties: 

92pii - 9ipi2 and 9ip22 - 92pi2 = 0. 

These kinds of observation will lead us to update Table [T] to give Table [H We have used the notation tt„ to 
denote projection onto an n-dimensional irreducible submodule. 

Now consider the algebraic Jacobian in genus two. This is a locus of quadratic identities in P^^. The 
curve is, as before, 

2/2 = Xqx^ + 6A5X^ + 15A4a-'' + 20X3X^ + 15\2X^ + 6Xix + Ao 

Equivariant inhomogeneous coordinates on P^^ are symmetric functions in a general divisor ((xi, yi), {x2, 2/2)) 
of degree two on the curve having poles of order two or more on the divisor. There are fifteen coordinates 
which fall into irreducible 5*^2 (C) modules of dimensions one to five: 

15 = 52 © 43 ® 34 ® 25 © Ig 

It is a crucial observation that the divisor pole grading (denoted by the subscripts) is simply related to the 
module dimension, i.e. that the dimension plus the degree at infinity is always equal to 7. 
The explicit coordinates on Jac{X) are: 



52 = 

43 = 

34 = 

25 = 

l6 - / 



1 2(a;i+a:;2) x\+4:XiX2+x\ 2a;ia;2(a;i + X2) x\x\ 
yi - y2 X2y\ - 0:12/2 xlyi - xly2 x^yi - xly2 

(53 ' 53 ' S3 ' §3 

I I{xi + X2) 1x1X2' 



(5' 5 ' 6 

yilyxi +y2l,X2 X2yilyxi +Xiy2lx2 
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where 6 ~ xi~X2 and S^I is equal to the equivariant polar form of the curve, discussed earlier. The vanishing 
of this polar form defines a cubic curve with order three contact at a point on the sextic curve above. 

These coordinates are equivariant modifications of coordinates introduced in [15j and the Jacobian is the 
locus of a set of 72 identities quadratic therein. We seek such identities by tensoring up the coordinate mod- 
ules, decomposing into irredueibles, identifying the divisor pole grading and then balancing both singularities 
and dimension. For example the projection 

7r5[34 34 -52 16] 

onto a 5-dimensional component is, as is easily verified, vanishing. This represents 5 quadratic identities. 
The other 67 can be equally straightforwardly described and classified. [5]. 

The decompositions of tensor products of standard irreducible modules of dimension d and pole order 
i on the divisor satisfy relations of the following type 

I n — m I — 1 

lis Ulr G Pp(n,s;7Ti,j-,p) • 

p—\n—m\-\-l 

The function p{n, s\m,r,p) is the pole order of the irreducible component of dimension p and it is an 
interesting problem to understand its structure. 

These ideas also apply to tensoring of piii2...i„i which is why we could rewrite all of Table [T] in terms of 
modules in Table [51 By a simple counting argument one sees that 

dimr(m + 1) — dimr(TO) = dimy(m + 2) + dimT^(TO — 1), 
i.e. (m + 1)^ - = {m + 2) + (m - 1). 

The Hirota derivative, P, has been introduced earlier as 

and it can be characterised in terms of 5'L2-modules by the fohowing equivariance property: 

V(g)V — > V(g)V ^ V 

4- J' 
V(g)V — > V(S)V — > V 



In this (commuting) diagram the downward vertical arrow signifies the 5'L2-action. Extensions of this 
equivariant Hirota map to higher order tensor products exist and our current work seeks to understand basis 
generation (in "appropriate" gradings) by studying equivariant resolutions of S'L2-modules of functions on 
the Jacobian. 
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